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ABSTRACT 

• We re-examine the existence and extent of the planar structure in the local galaxy 

t^j- \ density field, the so-called Supergalactic Plane (SGP). This structure is studied here 

ff") , in three dimensions using both the new Optical Redshift Survey (ORS) and the IRAS 

Q\ • 1.2 Jy redshift survey. The density contrast in a slab of thickness of 2Qh~ 1 Mpc and 

diameter of 80 ft. -1 Mpc aligned with the standard de Vaucouleurs' Supergalactic co- 
ordinates, is 5 S gp ~ 0.5 for both ORS and IRAS. The structure of the SGP is not 
well described by a homogeneous ellipsoid, although it does appear to be a flattened 
^ <— | structure, which we quantify by calculating the moment of inertia tensor of the density 

O field. The directions of the principal axes vary with radius, but the minor axis remains 

within Z ~ 30° of the standard SGP Z-axis, out to a radius of 80 /i" 1 Mpc, for both 
ORS and IRAS. However, the structure changes shape with radius, varying between a 
flattened pancake and a dumbbell, the latter at a radius of ~ 50 ft." 1 Mpc, where the 
| Great Attractor and Perseus-Pisces superclusters dominate the distribution. This calls 

to question the connectivity of the 'plane' beyond ~ 40 ft. -1 Mpc. The configuration 
found here can be viewed as part of a web of filaments and sheets, rather than as 
an isolated pancake-like structure. An optimal minimum variance reconstruction of 
the density field using Wiener filtering which corrects for both redshift distortion and 
shot noise, yields a similar misalignment angle and behaviour of axes. The background- 
independent statistic of axes proposed here can be best used for testing cosmological 
models by comparison with TV-body simulations. 
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1 INTRODUCTION 

The major planar structure in the local universe, the Su- 
pergalactic Plane (SGP), was recognized by de Vaucouleurs 
(1953, 1956, 1958, 1975a,b) using the Shapley-Ames cat- 
alogue, following an earlier analysis of radial velocities of 
nearby galaxies which suggested a differential rotation of 
the 'metagalaxy' by Rubin (1951). This remarkable feature 
in the distribution of nebulae was in fact already noticed by 
William Herschel more than 200 years earlier (for historical 
reviews see Flin 1986, Rubin 1989). 

When referring to the SGP in this paper we mean the 
planar structure in the galaxy distribution. This should be 



distinguished from the formal definition of Supergalactic 
coordinates, de Vaucouleurs et al. (1976, 1991) defined 
a spherical coordinate system L, B in which the equator 
roughly lies along the SGP (as identified at the time), with 
the North Pole (B — 90°) in the direction of Galactic coordi- 
nates (I = 47.37°; b = +6.32°). The position L = 0°,B = 0° 
is at (I = 137.37°, b — 0°), which is one of the two regions 
where the SGP is crossed by the Galactic Plane. The Virgo 
cluster is at (L = 104° ;B = -2°). Traditionally, the Virgo 
cluster was regarded as the centre of an overdense region 
called the 'Local Supercluster' (e.g. Davis & Huchra 1982, 
Hoffman & Salpeter 1982, Tully & Shaya 1984, Lilje, Yahil 
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& Jones 1986). However, some of the much larger super- 
clusters that are seen in recent redshift surveys, such as the 
Great Attractor and Perseus-Pisces, are near the SGP, and 
are possibly connected with the Local Supercluster, in the 
sense of being simply connected by an isodensity contour 
above the mean density (cf., Strauss et al. 1992; Strauss & 
Willick 1995; Santiago et al. 1995). 

Although the SGP is clearly visible in whole-sky galaxy 
catalogues (e.g. Lynden-Bell & Lahav 1988; Lynden-Bell et 
al. 1988; Raychaudhury 1989; see also the references above), 
the extent of planar structure in the galaxy distribution in 
3 dimensions has seen little quantitative examination in re- 
cent years. Tully (1986, 1987) claimed that the flattened 
distribution of clusters extends across a diameter of ~ 0.1c 
with axial ratios of 4:2:1. Shaver & Pierre (1989) found that 
radio galaxies are more strongly concentrated to the SGP 
than are optical galaxies, and that the SGP as represented 
by radio galaxies extends out to redshift z ~ 0.02. Stanev 
et al. (1995) claimed that energetic cosmic rays arrive pref- 
erentially from the direction of the SGP, where potential 
sources (e.g. radio galaxies) are concentrated, but this re- 
sult was criticized by Waxman, Fisher & Piran (1997). Di 
Nella & Paturel (1995) revisited the SGP using a compi- 
lation of nearly 5700 galaxies larger than 1.6 arcmin, and 
found indeed that galaxies were concentrated close to the 
Supergalactic Plane. Loan (1997) and Baleisis et al. (1998) 
searched for the SGP in projection in the 87GB (north) and 
PMN (south) radio surveys, and found signatures at the 1 
and 3 a levels, respectively. 

The existence of a pancake-like structure has important 
theoretical implications. Gravity amplifies deviations from 
sphericity, i.e. an initial oblate homogeneous ellipsoid evolves 
into a disk, and an initial prolate structure ends up as a spin- 
dle (Lin, Mestel & Shu 1965). White & Silk (1979) modeled 
the Local Supercluster as a homogeneous ellipsoid in ex- 
panding universe, and considered implications for cosmolog- 
ical models and initial conditions. Although such a simple 
model is insightful, the structure of the SGP is far more com- 
plicated than a homogeneous ellipsoid, as we show below. In 
a more realistic cosmological scenario, where the primor- 
dial density field is Gaussian, non-spherical shapes are more 
abundant than spherical ones (Doroshkevich 1970, Bardeen 
et al. 1986), and hence are expected to appear at the present 
epoch as even more non-spherical. In the context of the 'top- 
down' Hot Dark Matter scenario, Zel'dovich (1970) showed 
that pancake-like structures are the natural outcome of grav- 
itational instability in the quasi-linear regime. Further stud- 
ies have indicated that a web of filaments could also form in 
hierarchical Cold Dark Matter (bottom-up) scenarios (e.g. 
Bond, Kofman & Pogosyan 1996), but the shapes look quite 
different in different scenarios, as a result of both initial con- 
ditions and the cosmic time which allows the perturbations 
to grow. Hence, quantitative measures of the SGP and other 
observed filamentary structures and superclusters (Bahcall 
1988) could be very important in distinguishing between 
models, e.g. by applying the shape statistic to both data 
and to iV-body simulations. 

Here we study the properties of the SGP using the Op- 
tical Redshift Survey (ORS, Santiago et al. 1995) which 
provides the most detailed and uniform optically-selected 
sample of the local galaxy density field to date. For com- 
parison, we also use the full-sky redshift survey of galax- 



ies selected from the database of the Infrared Astronomical 
Satellite (IRAS), complete to 1.2 Jy at 60^im (Fisher et al. 
1995a). In this paper we consider the approach of 'moment 
of inertia' (Mol) to quantify a planar structure (cf. Babul 
& Starkman 1992; Luo & Vishniac 1995; Dave et al. 1997; 
Sathyaprakash et al. 1998). 

The outline of the paper is as follows. Section 2 de- 
scribes the ORS and IRAS samples, while in Section 3 
we discuss the visual impression and preliminary analysis 
of the SGP. Section 4 describes the formalism of the Mol 
we also develop related statistics which are independent of 
the background determination. Section 5 gives the results 
of Mol when applied to the samples, and Section 6 presents 
Wiener reconstruction of the Mol for IRAS. Section 7 gives 
interpretation of the results using mock realizations, and 
conclusions and future work are discussed in Section 8. 



2 THE ORS AND IRAS SAMPLES 

The Optical Redshift Survey (ORS, Santiago et al. 1995) 
covers the sky at Galactic latitude |6| > 20°. The survey was 
drawn from the UGC (Nilson 1973), ESO (Lauberts 1982), 
and ESGC (Corwin & Skiff 1995) galaxy catalogues, and it 
contains two subsets: one complete to a blue photographic 
magnitude of 14.5, and the other complete to a blue major 
diameter of 1.9 arcmin. The entire sample consists of 8457 
galaxies, with redshifts available for 98% of them; w 1300 
new redshifts were measured to complete the survey. The 
high number density of galaxies in ORS makes it ideal for 
cosmographical studies of the local universe. 

As ORS only covers \b\ > 20°, we filled in the Zone of 
Avoidance (ZOA) at \b\ < 20° with galaxies from the IRAS 
1.2 Jy survey (Fisher et al. 1995a). The Zone of Avoid- 
ance in IRAS, \b\ < 5°, was filled by interpolation based on 
the observed galaxy distribution below and above the ZOA 
(Yahil et al. 1991). While in principle it is possible to in- 
terpolate for the ORS ZOA, e.g. by hand or using a Wiener 
reconstruction (Lahav et al. 1994), we prefer to include real 
structure as probed by IRAS even at the price of mixing two 
different data sets. Hereafter when we refer to the ORS sam- 
ple we mean that supplemented by the IRAS 1.2 Jy galaxies 
at \b\ < 20°. 

Due to Galactic extinction and the diversity of cata- 
logues used, the selection function of ORS depends on both 
distance and direction (Santiago et al. 1996). In order to 
account for these selection effects, a weight is usually associ- 
ated with each galaxy. For a uniform survey like IRAS (see 
below) the selection function <f> depends only on the distance 
to a galaxy, |x|, and hence the weight is simply 

W9al = WW)' (1) 

where the mean number density of galaxies is estimated by 

<»> = p:£lMM) (2) 

gal 

(see Davis & Huchra 1982 for an alternative minimum vari- 
ance estimator of the mean density). 

For ORS, we follow Santiago et al. (1996) and Hermit 
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et al. (1996) and weight each galaxy by: 



Wgal 



1 (ri/ fj,,iras 

(n) M 0(x) (n)tot,ira 



(3) 



where fi labels the catalogue from which the galaxy in ques- 
tion is drawn; each of the four catalogues making up the 
ORS has its own mean density (n) M . Here (n) Mi i r . as is the 
mean number density of IRAS galaxies in the volume corre- 
sponding to the ^i-th catalogue in ORS. This brings each of 
the subsamples, with their separate selection criteria, to a 
uniform weighting (Santiago et al. 1996). With the weight- 
ing of equation (1) or (3), we can thus calculate the galaxy 
density fluctuation field <5(x) = (p(x) — (p))/{p) throughout 
the volume surveyed. 

For comparison we also apply our analysis to the com- 
plete IRAS 1.2 Jy sample (Fisher et al. 1995a; Strauss et al. 
1992), which includes 5313 galaxies over 87.6 % of the sky. 
The ZOA in IRAS was interpolated according to the Yahil 
et al. (1991) procedure mentioned above. Hereafter we re- 
fer to this sample as the IRAS sample. We note that the 
ORS and IRAS samples are not really independent; roughly 
60% of the IRAS 1.2 Jy galaxies out to 80 h' 1 Mpc are also 
catalogued in the ORS magnitude limited sample. 

We work purely in Local Group redshift space, except 
in section 6 on Wiener reconstruction. 




-30 30 -30 30 
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3 A BRIEF TOUR OF THE SGP 
3.1 Visual impression 

Historically, the SGP was identified in projected 2-D maps. 
We used the IRAS sample (which is homogeneous over the 
sky) to revisit the 2-D appearance of the SGP by fitting the 
data to a great circle. For a projected distribution of sources 
the great circle along which the density of galaxies is en- 
hanced can be found by calculating the covariance matrix: 
hj — ^2 ga[ XiXj where the Cartesian axes Xi (i = 1, 2, 3) are 
defined for a unit sphere. The matrix Uj can be diagonalised 
(cf. Section 4.2) and the smallest eigen-vector indicates the 
direction of the normal to the great circle. We considered the 
projected IRAS galaxy distribution out to 40 hT 1 Mpc and 
80 h -1 Mpc, and found that the great circle is aligned with 
the standard de Vaucouleurs' SGP to within 19° and 7° re- 
spectively. Can we see the translation of the great circle to a 
plane in 3-dimensions? Figure |l] shows the raw distributions, 
uncorrected for selection effects of ORS and IRAS galaxies 
in a sphere of radius of 40 ft -1 Mpc projected on the stan- 
dard (SGX-SGZ), (SGX-SGY) and (SGY-SGZ) planes. 
The SGP is particularly visible edge-on in the (SGX-SGZ) 
and (SGY-SGZ) projections. We note that the structures 
seen in ORS and IRAS are quite similar, but the ORS map 
is much denser, and clusters are more prominent. Maps of 
the smoothed density field, after applying the appropriate 
weights, are shown elsewhere for the IRAS survey (Saunders 
et al. 1991; Strauss et al. 1992; Fisher et al. 1995b; Strauss 
& Willick 1995; Webster, Lahav & Fisher 1997) and for ORS 
(Santiago et al. 1995, 1996; Baker et al. 1998). In Figures | 
and [] we show histograms of the galaxy density field 8, cor- 
rected for selection effects, as a function of SGZ, averaged 
over SGX and SGY, within spheres of ever-larger radii, as 
indicated in each panel. The SGP is seen as an enhancement 



Figure 1. The distributions of ORS and IRAS galaxies in a 
sphere of radius of 40 ft -1 Mpc, projected on the standard (SGX- 
SGZ), (SGX-SGY) and (SGY-SGZ) planes. The SGP is visible 
edge-on as the linear feature at SGZ = in the (SGX-SGZ) and 
(SGY-SGZ) projections. 



at SGZ = in these plots to roughly R = 60 h' 1 Mpc, 
but is not apparent on larger scales. 



3.2 The level of overdensity 

Our visual impression is that there is indeed a flattened 
structure in the galaxy distribution, aligned along SGZ = 0, 
extending to appreciable redshifts. To quantify this we cal- 
culate the overdensity S sgP = n sgp /{n) — 1 as a function 
of R m ax, where Ti S gp is evaluated within a circular cylinder 
centered on the Local Group with radius R m ax and height 
along SGZ of 20 ft" 1 Mpc, and the mean density (n) is cal- 
culated within a sphere of radius R ma x (including the slab). 
This overdensity reaches a maximum at Rmax = 40 h~ Mpc 
(for which the volume of the slab is ~ (46 ft -1 Mpc) 3 ), 
<W40) = 0.48 ± 0.06 and 0.46 ± 0.09 for ORS and IRAS, 
respectively. At R max = 60 h" 1 Mpc, 8 agp (60) = 0.36 ± 0.06 
and 0.26 ± 0.08 for ORS and IRAS, respectively. For refer- 
ence, the rms fluctuations in cubes of volume (40 hT 1 Mpc) 3 
are <5 rms ~ 0.37 and ~ 0.30 for optical (Stromlo/APM) and 
IRAS (1.2 Jy) surveys, respectively (Efstathiou 1993). Hence 
the fluctuation in galaxy counts in a slab aligned with the 
standard de Vaucouleurs' SGP is no more than ~ 1.5 o per- 
turbation on scale of ~ 40 /i -1 Mpc. This implies that the 
SGP is only a modest perturbation on these large scales. 
The overdensity as a function of the Z component of the 
standard SGP is depicted in Figures 2 and 3. 
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Figure 2. Histograms of the ORS density field (corrected for 
selection effects) as a function of SGZ, averaged over SGX and 
SGY, within spheres of radius as indicated in each panel. The 
error bars are Poissonian. 



3 L I ' 1 ' 1 I ' 
2 I R ™= 10 A "'Mpc 




-10 -5 
> M i | i i i I i 
: R m „=30 i-'Mpc 



10 -20 -10 10 SO 
R m „J40 h^Mpc ' 




-40 -20 20 40 



h K**"" e0 *~'llpo 4 



-40 -20 20 40 



1 h |i 

5 : T R ...t= 70 h Z M P° 

~-l 

-0.5 
I 




0.5 P 

~ 

0.5 i- 
-1 



R,.,,, jHO /HMp. 



-50 



50 -50 

SGZ (h-' Mpc) 



Figure 3. Histograms of the IRAS density field (corrected for 
selection effects) as a function of SGZ, averaged over SGX and 
SGY, within spheres of radius as indicated in each panel. The 
error bars are Poissonian. 



3.3 The centre of the overdensity 

In order to quantify the local overdensity, we need to calcu- 
late the centre of 'mass' of the density field, within a sphere 
of radius Rmax centred on us, for both the ORS and IRAS 
samples (cf., Juszkiewicz et al. 1990). The centre of mass is 
given by summing over the galaxies (assuming they all have 
equal mass): 



- - 1 V" 

&i — -yr / ^ IMgal %i j 
gal 



(4) 
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(b) IRAS 



R m .„ (mpc h-«) 

Figure 4. The variation of the mean position of the centre of 
mass of the galaxy distribution within a sphere of radius R m ax 
centred on us, for both the ORS and IRAS samples. The curves 
show the Supergalactic x, y, z of the centre of mass, and its radial 
distance from the origin. 



where i = 1,2,3, V is the total volume of the sphere of 
radius Rmax and w ga i is the weight per galaxy described in 
82. The shot noise error in each axis is 



1 

V 1 



E2 2 
Wgal Xi = 



gal 



1 

3V^ 



gal 



2 

'gal 



(5) 



where the last equality is valid for an isotropic distribu- 
tion. Figure ^ shows the coordinates of the centre of the 
galaxy distribution as a function of Rmax ■ The centre moves 
away from the origin (the Local Group) by no more than 
8/i _1 Mpc for ORS and 6/i _1 Mpc for IRAS. This is par- 
tially due to the 'tug of war' between the Great Attrac- 
tor and Perseus-Pisces, which largely balance each other 
out. Indeed, the centre moves back towards the origin for 
R > 40 /i" 1 Mpc. 



4 PRINCIPAL COMPONENTS OF THE 
INERTIA TENSOR 

Both visual impression of and a formal y 2 fit suggest that the 
SGP cannot be fitted by a homogeneous ellipsoid model. We 
therefore seek a more objective measure of deviation from 
sphericity, via the Moment of Inertia tensor (Mol) . 



4.1 Estimation of the Moment of Inertia 

We wish to detect a high density planar structure (e.g. a 
slab or an ellipsoid) embedded in a uniform sphere of radius 
Rmax- One approach (cf. Babul & Starkman 1992, Luo & 
Vishniac 1995, Dave et al. 1997) is to construct the Mol for 
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the fluctuation in the density field: 



Cij Cij ' 



Cij — 



1 



(n)V 



i(x.)-n bg ](xi-Xi)(xj-Xj) dV, 

(6) 

where n(x) is the galaxy number density at position x, (n) 
is the mean number density and Xi,Xj = 1,2,3) are 
Cartesian components of x. Note that we allow the centre 
x to move (see Figure 4), although below, we will end up 
keeping it fixed. The integration is over the volume of the 
sphere. We define the fluctuations in the density field relative 
to the background density n;, ff in the absence of the slab 
(which differs from the mean density (n) including the slab). 

For a uniform distribution with density (n) the covari- 
ance matrix is analytic and isotropic: 



/=< _ fbg 



1 pK fbg 

3 ij V 



p2 



r (4ivr dr) 



&ij fi 



bg 



(7) 



space (e.g. Murtagh & Heck 1987 and references therein). By 
identifying the linear combination of input parameters with 
maximum variance, PCA determines the Principal Compo- 
nents that can be most effectively used to characterize the 
input data. In our case, in searching for a plane, we wish to 
see if the PCA finds one axis to be much shorter than the 
other two axes. 



4.3 Background-independent statistic of axes 

Our aim is to quantify deviations from spherical struc- 
ture. The difficulty in doing so is due to two problems: (i) 
the background density as modeled above with fbg is ill- 
determined, and (ii) the shot noise due to the finite number 
of galaxies can be quite large. 

We will discuss the shot noise problem in §6 below. To 
overcome the problem of the unknown background, we can 
construct the following quantities which are independent of 

fbg- 



where 



(12) 



r "69 

Ib3 ~ <n> ' 



(8) 



x is of course zero, and is the Kronecker delta function. 

For a discrete density field, such as provided by ORS or 
IRAS , we can calculate the covariance elements by summing 
over the galaxies: 



Cij — Cij 



Cij — — W ga lXiXj 
gal 



Rn 



ij Jbg - 



(9) 



where we have set x = 0, given the fact that the centre 
of mass stays so close to the origin (Figure 4). An analytic 
estimator for Poisson shot-noise per diagonal component of 
the covariance matrix (assuming no errors in fb g ) is: 



(31/)' 



E2 4 
Wgal T 



(10) 



gal 



Eq. (9) requires knowledge of fbg in order to correct for 
the background density. However, as the boundaries of the 
SGP are hard to define a priori and the background will 
not be uniform, it is difficult to estimate a meaningful fb g . 
In §4.3 we suggest a statistic which is independent of this 
parameter. 



4.2 Principal axes 

The next step in our analysis is to diagonalize the Mol and 
find the eigenvalues X a and associated eigenvectors u a (a = 
1,2,3): 



(ii) 



The standard deviation ('1-<t') along each of the three axes 
is given by which we label hereafter as a, b, c. Note that 
since the background contribution (the last term in eq. 9) 
is isotropic, it only affects the eigenvalues, but not the di- 
rections of the eigenvectors. This procedure is essentially 
the Principal Component Analysis (PCA) - a well known 
statistical tool for reducing the dimensionality of parameter 



2 I 2 2i 

q — \a — c | 



and 



(13) 



(14) 



For a perfect sphere, p — q = s — 0. For a homogeneous 
oblate ellipsoid with semi-major axes A = B > C, one ex- 
pects p = for all R m ax, but for q and s to increase out 
to Rmax ~ A and then to decline. Note that the (p, q, s) 
here are in the diagonalized PCA frame, where the axes are 
ordered by size (a > b > c). Alternatively, we can calculate 
(p, q, s) along fixed (x, y, z) axes (e.g. in de Vaucouleurs' 
system) by replacing a 2 ,b 2 ,c 2 in eqs. (12-14) by the non- 
ordered variances in that fixed coordinate system (i.e. the 
diagonal elements in the Mol matrix in that (x, y, z) coor- 
dinate system). 

In the limit that the errors are isotropic the Poisson 
error bars are: 



= o 2 c /p . 



(15) 



and similarly for q and s. We will determine the character- 
istic shape of the galaxy density field, from the observed 
p, q and s and their error bars as a function of Rmax ■ We 
note that other transformations of (a, b, c) are possible, such 
as the (SI, 52, S3) statistic proposed by Babul & Starkman 
(1992). 



4.4 Mock realizations 

To get an insight into our (p, q, s) statistic we applied it 
to mock realizations of ellipsoids placed in a uniform back- 
ground. The IRAS selection function was applied to the 
mock samples, and the results were averaged over 100 re- 
alizations. Figure [| shows the quantities p, q, and s derived 
at fixed axes for mock ellipsoids (oblate, prolate and triaxial) 
with dimensions indicated in the panels and overdensity of 
5 — 0.4. When the axes were allowed to be chosen by PCA, 
the axes agreed with the correct orientation of the ellipsoids 
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Figure 5. The quantities p (crosses), q (open squares) and s 
(filled circles) given by eqs. (12-14) vs. Rmax for mock ellipsoids 
(oblate, prolate and triaxial) with dimensions (in h^ 1 Mpc) indi- 
cated in the panels, all ellipsoids with density contrast of S = 0.4. 
The mock samples have an IRAS selection function, and the re- 
sults are averaged over 100 realizations at fixed axes. 



to within 7° (rms over 100 realizations). For a single real- 
ization the curves are much more noisy. Moreover, if PCA 
defines the axes, it attempts to maximize the differences be- 
tween the axes, and the results for p, q, and s are biased by 
noise. For example, in the case of a perfectly oblate struc- 
ture in the presence of noise, the quantity s increases with 
radius Rmax, where it is identically zero in the noise-free 
case. Fortunately, we will see in the following section that 
the behaviour of (p, q, s) for the real data is quite robust. We 
will continue these experiments with the mock realizations 
in §7, where we consider more complex geometries. 



particular in ORS), more in line with the visual impression 
of a pancake-like structure in Figure 1. However, the dra- 
matic changes with Rmax (in particular at 50 h^ 1 Mpc) call 
to question the notion of a single coherent feature within the 
boundaries of our samples (see §7). Based on Figure 6, one 
may argue that the 'plane' terminates at ~ 40 K~ x Mpc. 

We now turn to the PCA approach of diagonalizing 
the Mol derived at each radius Rmax- In this case the axes 
change direction, and might even 'swap'. 'Swapping' means 
that an axis which is the largest at a given Rmax may be- 
come say the second or third largest in a different Rmax- 
Hence we may see a discontinuity in the variation of the 
'largest' axis with R m ax- 

The variation at large Rmax is more difficult to inter- 
pret, but we see that q > p > s, implying a triaxial shape, 
and p and q increase with R ma x- In contrast, the toy models 
(Figure 5) show that two of [p, q, s] drop with Rmax beyond 
the radius which encloses the ellipsoid (80 ft. -1 Mpc in the 
toy models). But in making this comparison with Figure 5 
we should keep in mind that at large Rmax the increase of 
the observed quantities p, q or s might be due to shot-noise. 
It is also interesting to note the similarity between ORS and 
IRAS, indicating that on the very large scales the density 
fields in the two samples are quite similar. 

Table 1 lists a, b, c, their Poisson errors, and the direc- 
tion of the PCA Z component for ORS (for direct sum- 
mation) and IRAS (for both direct summation and Wiener 
reconstruction described in the next Section). Table 1 also 
lists the values in differential shells, showing trends similar 
to those seen in the cumulative plots. 

Fig 9 shows the variation of direction of the PCA Z 
component with Rmax- We see that although the angle 
9 Z = 90° — B z between the PCA minor axis and de Vau- 
couleurs' pole is usually less than 30°, there is a significant 
variation on the sky (relative to the errors of ~ 10° due 
to shot noise). Particularly noticeable is the big jump at 
Rmax = 50 hT 1 Mpc, in accord with the behaviour in Figures 
6,7 and 8. This is again due to the dumbbell configuration 
of the Great Attractor and the Perseus-Pisces superclusters, 
as we shall see in §7. Table 1 also gives B z , the complement 
of the misalignment angle 8 Z , for 3 shells, each 20 hT x Mpc 
wide. We see that indeed typically 8 Z ~ 30°, even at the 
(60-80) h^ 1 Mpc shell, which is beyond the Great Attrac- 
tor and Perseus-Pisces. The probability that two random 
vectors are separated by an angle less than 6 Z ~ 30° by 
chance is P{< 6 Z ) = [1 - cos(0 z )], i.e. ~ 13%. 



5 APPLICATION OF MOI TO ORS AND IRAS 

We now apply the Mol approach to the ORS and IRAS sam- 
ples. We begin by considering p, q, and s along fixed axes. 
We will find below that the PCA direction at small Rmax is 
indeed very close to that of de Vaucouleurs' system, so wc 
simply choose the standard de Vaucouleurs' SGP axes. Fig- 
ure 6 shows p, q and s (derived here from the variances along 
the standard SGP axes without ordering them, as explained 
in section 4.3) as a function of Rmax - The most notable fea- 
ture is the dip in s at about R max = 50 /i -1 Mpc in both 
samples. This is due to the dumbbell-like configuration of 
Perseus-Pisces and the Great Attractor on opposite sides 
of the sky, as we'll see in §7. But at R max = 80 h~ x Mpc, 
the structure is pancake-like: p tends to small values (in 



6 WIENER RECONSTRUCTION 

In addition to the problem of shot-noise, the density field 
measured from redshift surveys also suffers from redshift 
distortions due to peculiar velocities (see Hamilton 1998 for 
a review). One approach which deals with both problems 
is Wiener filtering (e.g. Lahav et al. 1994, Fisher et al. 
1995b). In particular, Fisher et al. (1995b) expanded the 
density field 8(r) in terms of spherical harmonics and Bessel 
functions. By assuming a prior model for the power spec- 
trum, the coefficients of the expansion can be corrected for 
the redshift distortion and for the shot-noise. This recon- 
struction is optimal in the minimum variance sense. In this 
approach, the smoothed density field approaches the mean 
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20 40 60 80 100 



R mH (Mpc h">) 

Figure 6. The quantities (p, q, s) derived from the variances 
along the de Vaucouleurs' axes vs. Rmax (cumulative) for the 
ORS and IRAS samples. The error bars are Poisson. 

density at large distances. This does not mean necessarily 
that the SGP itself disappears at large distances; it merely 
reflects our ignorance of the density field where the data are 
poor. 

This Wiener approach applied to the density field S also 
gives the optimal reconstruction for any property which is 
linear in S. In particular, if we seek the optimal reconstruc- 
tion of the moment of inertia (eq. 7) it can be re- written (for 
x = 0) as 

^^stl^^-A'l?^ ' (16) 

where 

hj = I 8{Y)xiXjdV . (17) 

J R 

Webster, Lahav & Fisher (1997) give analytic expressions 
for Iij in terms of the reconstructed coefficients 5{^ nn , and 
show that only harmonic modes I = and I = 2 contribute. 

We apply this technique below to the reconstructed, 
real-space density field of IRAS galaxies, assuming as priors 
P = Q? 6 /b — 0.7, normalization erg = 0.7 and a CDM power 
spectrum with shape parameter F — 0.2. 

Figure 10 shows the p, q and s for IRAS before and 
after Wiener filtering in spherical harmonic presentation. 
Figure 11 shows the variation of the minor PC A axis with 
Rmax- There is good agreement with the results of direct 
summation shown in Figure 9, indicating that shot noise and 
redshift distortion do not introduce big systematic effects in 
our statistics. Unfortunately, the Wiener procedure cannot 




20 40 60 80 100 



R mH ("PC M 

Figure 7. The quantities (p, q, s) derived from PCA dimensions 
(a, 6, c) (eqs. 12-14) vs. Rmax (cumulative) for the ORS and IRAS 
samples. The error bars are Poisson. 




-1 -0.5 0.5 1 

SGX = cos(SGB) cos(SGL) 

Figure 8. ORS : The direction of the Z axis from the PCA 
analysis (using direct summation over the galaxies) relative to 
the de Vaucouleurs' pole. Values of Rmax arc indicated. 



easily be applied to the ORS sample, as it only covers |b| > 
20°. 
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-1 -0.5 0.5 1 

SGX = cos(SGB) cos(SGL) 

Figure 9. IRAS : The direction of the Z axis from the PCA 
analysis (using direct summation over the galaxies) relative to 
the de Vaucouleurs' pole. Values of Rmax are indicated. 




- IRAS^.^iene-r;- ttgrntonics - 

_ 1 i i i i '' i "'r'^;:r^r'V' i i i i 

-1 -0.5 0.5 1 

SGX = cos(SGB) cos(SGL) 

Figure 11. Wiener filtering of IRAS : The direction of the Z axis 
(as found by PCA) relative to the de Vaucouleurs' pole. Values 
of Rmax are indicated. 
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Figure 10. The quantities p, q, s vs. Rmax (cumulative) in PCA 
axes using spherical harmonic reconstruction of IRAS and an- 
alytic Mol. Top: corrected for rcdshift distortion and noise by 
Wiener filtering, with parameters /3 = 0.7, T = 0.2, erg = 0.7 (re- 
construction out to 200 h^ 1 Mpc). Bottom : using raw harmonics. 



7 INTERPRETATION OF THE MOI RESULTS 

The behaviour of the Mol axes (Figures 6-10) for both ORS 
and IRAS indicate that the SGP is a structure far more com- 
plicated than a homogeneous ellipsoid or a slab. In the local 
universe there is 'contamination' due to clusters and voids 
below and above the SGP. To get further insight to the struc- 
ture that generates the observed Mol, we generated mock 
realizations of the IRAS sample with more complicated ge- 
ometries than the simple ellipsoids used in §4.4, trying to 
mimic the significant overdensities seen in the real universe. 
In particular, we put in: (i) an ellipsoid with dimensions of 
70 x 70 x 10( ir 1 Mpc) 3 and overdensity 5 = 1 aligned with 
the Supergalactic plane, (ii) a sphere of radius 7.5 ft -1 Mpc 
and overdensity S = 2.5, representing a Virgo-like cluster 
at (SGX, SGY, SGZ) = (5, 10, -2) h' 1 Mpc, (iii) a sphere 
of radius 11 hT 1 Mpc and overdensity 5 = 3, representing 
a Perseus- Pisces- like supercluster at (SGX, SGY, SGZ) = 
(49, -23, 0) h' 1 Mpc, and (iv) a sphere of radius 10 h^ 1 Mpc 
and overdensity 5 = 3, representing a Great Attractor-likc 
structure at (SGX, SGY, SGZ) = (-32, 0, 0) ft" 1 Mpc. 

Figure 12 shows (p, q, s) averaged over 100 realizations 
in fixed and PCA coordinates. It is interesting to contrast 
this Figure with Figures 6, 7 and 10. For example we see 
that the dips in s at ~ 25 h -1 Mpc 50 hT 1 Mpc in the real 
data (Figure 6) are reproduced here (albeit with not quite 
the same amplitudes), especially with fixed axes. This toy 
example illustrates that the observed behaviour of p, q, and 
s can be accounted for by adding few major superclusters 
to a pancake-like structure. We note that the toy model in 
Figure 12 does not reproduce the observed growth of p, q 
and s on large scales, which might be due to shot noise. 

We note that recent simulations for a variety of models 
(e.g., Bond et al. 1996, Jenkins et al. 1998) show a web 
of filaments and sheets, with 'knots' at their intersection, 
or alternatively, 'arms stretching from clusters'. It may well 
be that Perseus-Pisces and the Great Attractor represent 
such knots. Hence the Mol 'shape-finder' statistic, when ap- 



© 0000 RAS, MNRAS 000, 000-000 



SGP parameters from PCA analysis 
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Table 1. 
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Comments on Table 1: R m i n , Umax, o,,b and c arc in h^ 1 
Mpc. The values of a, b and c arc for f og = 0. The Poisson error 
per axis A(a 2 ) is derived using eqs. (15-17). No errors are given 
in the Wiener case, as this is a smoothed field (although one can 
calculate the scatter around the mean). 

L z and B z (both in degrees) give the direction in standard de 
Vaucouleurs' coordinates of the Z component as found by PCA. 
The angle between this direction and de Vaucouleurs' SGP pole 
is simply 6 Z = 90° - B z . 

plied to many points in space, could be an efficient way of 
detecting such patterns and discriminating between models. 



8 CONCLUSIONS 

We have presented a Moment of Inertia approach to study 
the existence and extent of the planar structure in the local 
galaxy density field, the Supergalactic Plane (SGP). 

Our main conclusions from the present analysis are: 

(i) The SGP is not well-modeled as a homogeneous el- 
lipsoid. 

(ii) The Moments of Inertia analysis of the galaxy den- 
sity field shows that the normal to the plane is within 
6 Z ~ 30° of the standard SGP Z-axis, out to a radius of 
80 ft" 1 Mpc, for both ORS and IRAS. However, the normal 
changes direction with R ma x, in particular at ~ 50 h^ 1 Mpc, 

© 0000 RAS, MNRAS 000, 000-000 
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R m „ h-') 

Figure 12. The quantities p,q,s vs. Rmax (cumulative) in fixed 
and PCA axes for a mock realizations of an oblate ellipsoid plus 
three clusters, as described in the text. The upper panel uses 
fixed supergalactic coordinates, while the lower panel uses the 
axes found by the PCA analysis. Results were averaged over 100 
realizations. 



due to the presence of Perseus-Pisces and the Great At- 
tractor. The shape too varies with R m ax, and in fact at 
50 h' 1 Mpc the SGP looks more like a dumbbell than a pan- 
cake. Based on Figures 6, 7 and 10 one may argue that the 
'plane' terminates at ~ 40 ft -1 Mpc. However, we note that 
the axes the shell (60 — 80) hT x Mpc are aligned with the 
standard SGP to within 30°. 

(iii) The density contrast in the slab with Rmax = 
40 ft" 1 Mpc and thickness of 20 /i" 1 Mpc centered on the 
Local Group and aligned along the Supergalactic axes is 
S sgp ~ 0.5 for both ORS and IRAS. 

(iv) The SGP axes and density contrast are similar for 
both ORS and IRAS. 

(v) An optimal minimum variance reconstruction 
(Wiener filtering in spherical harmonics representation), 
which corrects IRAS for both redshift distortion and shot 
noise, yields similar misalignment of angle and axes. 

(vi) It may well be that Perseus-Pisces and the Great 
Attractor represent 'knots' in a web of filaments and sheets, 
as seen in recent Af-body simulations. This complicated 
structure calls to question the notion of a single coherent fea- 
ture on Gpc scales (as claimed e.g. by Tully 1987, 1987), and 
the possible connectivity between the Great Attractor and 
the Shapley Supercluster at ~ 140 ft" 1 Mpc (e.g. Scaramella 
et al. 1992). 

Regarding the wider cosmological implications, the 
present analysis is as an example of applying the Mol at only 
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one point, centred at us. It is of interest to compare the SGP 
to other structures, e.g. the Great Walls seen in the CfA and 
SSRS surveys. One can calculate moments of inertia at many 
points in space and contrast the statistical behaviour of the 
axes with iV-body simulations. Results from volume limited 
subsets of -/V-body simulations show indeed that the Mol 
statistic quantifies e.g. differences in filamentary structure 
between Cold and Hot dark matter models (e.g. Dave et al. 
1997, Webster 1998). Unfortunately, the flux limited ORS 
and IRAS catalogues are too shallow to apply the statistic 
in other points in space. Volume limited subsets from new 
large redshift surveys (e.g. 2dF and SDSS) would be suitable 
for such shape-finding algorithms. 

It also is possible to improve the above analysis by 
calculating the moment of inertia for a pre-selected re- 
gion above a certain density-threshold. It would also be 
interesting to repeat the analysis for different morpholog- 
ical types and weighting schemes; for example, by calcu- 
lating luminosity- weighted, rather than number weighted 
galaxy counts. On the statistical side, other more general 
shape-finders can be applied, e.g. by utilizing the Minkowski 
functions (e.g. Mecke, Buchert & Wagner 1994; Sahni, 
Sathyaprakash & Shandarin 1998). 
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